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FROZEN AND ALMOST FROZEN STRUCTURES IN THE 
COMPRESSIBLE ROTATING FLUID 

OLGA S. ROZANOVA 


Abstract. We study a possibility of existence of localized two-dimensional 
structures, both smooth and non-smooth, that can move without significant 
change of their shape in a leading stream of compressible barotropic fluid on 
a rotating plane. 


1. Bidimensional model of compressible fluid 

We consider the system of barotropic gas dynamics in 2D on the rotating plane, 

e(5iU + (u-v)u + /:u) + vp(p) = o, p = Cp\ c>o, (i) 

9*0 + div(£>U) = 0, (2) 

for density g, vector of velocity U and pressure P, t > 0, a; G Here £ = IL, 
P = ^ ^ Q ^1 is the Coriolis parameter assumed to be a positive constant, 
7 G (1, 2) is the adiabatic exponent. 

Under suitable boundary conditions system implies conservation of mass, 

momentum and total energy. 

Many models of ocean, atmosphere and plasma are approximately two dimen¬ 
sional. In particular, in [I] a procedure of averaging over the hight in a three- 
dimensional model of atmosphere consisting of compressible rotating polytropic 
gas was proposed (see also [5]). 

"y — 1 

Let us introduce a new variable tt = P . For 7 r(t, x) and U(t, x) we obtain the 
following system: 

dtV -f (U • V)U -L ILU + coVtt = 0, dtir + (Vtt • U) -f (7 - l) 7 r divU = 0, 

with Co = Then we change the coordinate system in such a way that the 

origin of the new system x = (xi,X 2 ) is located at a point X(t) = {Xi{t), X 2 {t)) 
(here and below we use the lowercase letters for x to denote the local coordinate 
system). Now U = u-|-V, where V(t) = {Vi{t),V 2 {t)) = {Xi{t), X 2 {t))- Thus, we 
obtain a new system 

9tu + (u-V)u-LV + /£(u-LV)+coV^ = 0, (3) 

dtTT + (Vtt • u) -|- (7 — 1) tt div u = 0. (4) 

Given a vector V, the trajectory can be found by integration from the system 
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In our previous papers m. m we used this approach to find a position of atmo¬ 
spherical vortex, associated with a tropical cyclone. 


2. Local and bearing fields separation 

We assume that the pressure field 7T{t,xi,X2) can be separated into two part as 
TT = TTo -l- TTi, where tti is somewhat stronger, however more uniform than ttq. We 
call TTo the local field and tti the bearing field. 

Then we obtain from ((S]), dH) 


[9tu 4- (u • V)u -f £u -I- Co Vtto] 4- 


V 4- £ V 4- Co Vtti 


= 0 , 


[dtTTo 4- (Vtto • u) 4- (7 - 1) TTo div u] 4- [dt-Ki 4- (Vtti • u) 4- (7 - 1) tti div u ] = 0. 
If we assume that we can find the couple (u, tto) from the system 

9tu 4-(u • V)u 4- /lu 4- Co Vtto = (5) 

dt'Ko + (Vtto • u ) -4 (7 - 1) TTo div u = 0, ( 6 ) 

with a certain function Q{t, x), then we get a linear equation for tti, 

dtTTi + (Vtti • u) 4- Tridivu = 0, (7) 

which can be solved for any initial condition 7 ri( 0 ,x). Further, ([5]) and ([3]) imply 
V(t) +£V(t) 4 - coV7ri(t,x) = -Q{t,x). 


Now we set 

Q =-Co [Vtti (t,x) - Vtti( t,0)]. (8) 

Thus, we associate the couple (u,7ro) with the local field and the couple (V,7ri) 
with the bearing field. As we can see, the couple (u,7ro) is independent of the 
bearing field ”up to the function Q.” If the solution (u, tto) to system ([5]), ([6]) is 
found, we can find (V,7ri) from linear equations. 

If Q = 0, then the bearing field does not influence on the local field and in this 
sense we will talk on a complete separation of he bearing and local fields. Evidently, 
this will be only if tti is linear with respect to the space variables. If |(5| < 5 for 
sufficiently small (5 > 0, we can talk about a ”S- approximate” separation of fields, 
Q plays a role of discrepancy. This discrepancy is a measure of separability of the 
local and bearing fields. 

The position of the center of the moving coordinate system can be found from 
the following equation: 


X(t) 4- £X(t) 4- CoV tti (t, x) 


= 0 . 

x=0 


(9) 


As follows from the computer modeling made in [4], even for real meteorological 
data the position of center of tropical cyclone found by means of equation ([9]) is 
quite accurate. 
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3. Steady nonhomogeneous incompressible flow 


We look for a solution of the local field with special properties, namely, a steady 
divergence free solution. If the discrepancy (5 = 0, then this solution can be 
considered as a ” frozen pattern” into a leading stream. If the discrepancy is small, 
we can talk only on an ” almost frozen pattern”, since the right-hand side in equation 
®, that depends on the properties of the bearing field, influences the solution. 

Thus, let us assume that 


u does not depend of t, u(0) = 0 and div u = 0. 
This means that there exists a stream function ^(xi,X 2 ) such that 

u = 

Equations (O and ([5]) result 


(V^o-Vi$) = 0, (10) 

V)Vi$-f ZLV_L$ + coV7ro = 0. (II) 

We take the inner product of (HU and and get 

^^1^2 - ($ 2 . 1 )^)) = (12) 

The solution of (fT^ have to satisfy the identity 


V X ((V_L$ • V)Vi$) = 0. 
Equations (HU) and m are equivalent to 

j($,iv$n = o, 

and 


( 13 ) 

(14) 


J($,A$)=0, (15) 

respectively, where J is the Jacobian. The stream function $ has to satisfy both 
da and da, and for smooth $ the function ttq can be restored up to a constant 
as 


1 

TTo =- 

Co 


(‘i> 


X2^XiX2 ^Xi^X2X2 




( ^X2^X\X\ H“ ^X\^X\X2 “1“ ^^ 3 : 2 ) dx2 

{^x2^xiX2 ^xi^X2X2) dx\ H" H" ^3:1^3:13:2) dx2 

There are two evident classes of solution to da and da: 

$ = ^{xf + xl) 


■ ( 16 ) 


and 

4> = $(xi), * = 1,2, 

with arbitrary smooth function of one variable $. In particular, $ can be compactly 
supported. 

The first case corresponds to a steady vortex. In the meteorological model this 
pattern can be associated with a tropical cyclone in the mature stage of develop¬ 
ment. 

The second case corresponds to a shear flow and can be associated with an 
atmospheric front. 
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Remark 3.1. J/tto = const, then the condition m is eliminated and the problem 
can be reduced to the solution of the Euler equations. Even in this case possible 
steady solutions can be very complicate 0. 0> H- 


Remark 3.2. Equations ^4^ and (03j imply the Dubreil-Jacotin equation [5] 


A$ + 


^o(^) 

2(7 - l) 7 ro($) 


|v$p 


Fm 


if we take into account that ttq = 7ro(<&). The last property follows from ilOi) . The 
arbitrary functions 7ro(<&) and F{^) give initial distribution of density and vorticity. 


3.1. Algorithm of solution, the smooth case. 

Theorem 3.1. Let <i>(a;i,a; 2 ) be a smooth solution to the system 

|V$p=C?($), (17) 


= (18) 

with a differentiable function G and integrable function R. Then $ solves 
and therefore it is a part of solution to the system m, (W- 

Proof. Equations (El) and (fTSl) means that |V$p and A$ and respectively, 
are functionally dependent, this results (IT4)) and El- Thus, <i> solves the system 
(fT0|. (fTT|) together with ttq found by El- Namely, taking into account El and 
El we obtain 





2i?i($) + 2l^)xi dxiP 


- 2Ri($) + 2Z$),,dx2 = 

^ (|V$p - 2i?i($) - 2/$) + const, 
2co 

where Ri = J R{r])dr]. 

<J>0 


□ 


3.2. Relation with the eikonal equation. 

Proposition 3.1. Assume that a function ^ satisfies in a domain LI the standard 
eikonal equation 

|Vf|2 = l. (19) 

Then any differentiable monotone function $ = iF{^) satisfies equation with 


Proof. Proof is a direct computation. 


□ 
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4. Construction of localized frozen patterns 


4.1. The smooth case. Let H be a compact domain in with smooth boundary. 
We assume that there exist a couple of functions G($), i?(4>)) such that the stream 
function $ satisfies equations (ITTl) . (ITSl) for x S in the classical sense and 

$|aa = 0, V$|an = 0. (20) 

Proposition 4.1. Let d* G be a solution to ( [17[ 1, il8\) for x € Q with 

boundary conditions li2(J\) . Then the solution to & in the domain \ in the 
local coordinate system do not depend of the local field. 

Proof. Condition (l20ll means u|an = 0, therefore the velocity field can be extended 
smoothly to the whole plane R^ as zero. Then the solution to (O keeps its initial 
value for X ^ fl. □ □ 

Thus, the solution to problems (fT71) . (IT^ . (EUl) gives a ’’frozen pattern” inside 
the domain Ll. 

Let us show that the class of solutions satisfying the conditions of Proposition 
EH is not empty. The simplest situation is where is a disc of radius r. The 
solution to (EU with boundary value 0 = r is ^ = \/x'f + the differentiability 
fails only at the origin. Nevertheless, the solution to (l20)l based on (EU can be 
smooth everywhere in if we take F{f^) = . Further, we can take F = A(s), 

s = with any smooth monotone on [0,1] function A such that A(r^) = 0, 
= 0, i = 1,2. Further, A$(0 = = 

This example gives a variety of axisymmetric vortex structures. 

Remark 4.1. As follows from [^, Sec.2.3.3, these structures are nonlinearly stable 
with respect to smooth perturbations keeping zero boundary conditions. 

4.2. Non-smooth case. Now we consider non-classical solution to the system 
EH), El, allowing to construct the "frozen patterns” containing discontinuities. 


ds^ 


d^X 

ds'^ 


4.2.1. Generalized solution to the eikonal equation {1 7[ 1. It is well known that the 
problem of constructing a nonlocal theory of the Cauchy problem for nonlinear 
Hamilton-Jacobi equation (in particular, for El and El) inevitably leads to the 
necessity of introducing a generalized solution. A natural extension of the notion 
of the solution in the sense of ’’almost everywhere”, i.e. a locally Lipschitz contin¬ 
uous function satisfying the equation everywhere in the considered domain except 
possibly at the points of a set of zero measure. Generally speaking, this solution is 
not unique. In El it was introduced the following notion of generalized solution, 
which we formulate applying to our case. 

Let Br{y) = {x G R^ : \x — y\ < r}. We denote by Cipi^^^D) the totality of 
functions /(x) defined on a set D G R^ and satisfying the Lipschitz condition on 
the subset Br{y) fl D. It is known El] that /(x) has at almost every interior point 
ot D a differential and hence a gradient Vu. 

Further, let fl is a bounded domain in R^ and / G Cipig^{Ct). We say that /(x) 
belongs to the stability class E{Vi) if the following inequality is satisfied for every 
X, X -I- Ax, X — Ax G Bs{y) C B 25 C 12, (Ax 0): 


A^/ ^ /(x -b Ax) - 2/(x) -b /(x 

|AxP |AxP 


Ax) 


> —C{y, 6) = const. 
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Definition 4.1. A function $ G Cipig^{VL) fl E{Q?) is called a generalized solution 
of the Cauchy-Dirichlet problem 

|V$p=G($), eiao=0, (21) 

where G(<i>) is a smooth and almost everywhere positive function, if it satisfies 
equation almost everywhere in 12 and takes the boundary values. 

4.2.2. Generalized solution to the nonlinear Poisson equation m- We give the 
definition of the generalized solution following [12] . 

Definition 4.2. 4> G is called the generalized solution to the boundary 

problem 

A$ = cc G f2, ^ = 4>, X G dfl, 

with R G 1^2 (12), 4> G if 

4> - (/) G i?o(12) 

and 

J (V<f>,Vt!)dx = — J R{^)vdx 

n n 

for all V G Hglfd). 


4.3. Example of discontinuous frozen pattern. Now we construct a discon¬ 
tinuous solution to the system (flTll . (IT^ in the sense of Definitions 14.11 and 14.21 
satisfying zero boundary conditions. 

We consider the domain 12 such that (xi, X 2 ) G 12 if xf -|- x| < 1 and Xi < i. 

It can be readily checked that the function 


^(X1,X2) 


2xi, —\/3xi < X 2 < V^xi, 

i/xf + x|, otherwise, 


solves the eikonal equation (IT^ in 12 with boundary value ^ = 1. Let us take again 
as F = A(s), s = l^^l, any monotone smooth function on [0,1] such that A(l) = 0, 

414 =0, i = l,2. 

ds^ 9 ds^ r, ’ ’ 

s—r‘^ s—0 _ 

Thus, $ G £ipig^{fl) solves (fT71) and (fT51) in the sense of Definition 14.21 with 
R{^) = $"(^) -b — y(Di)), where x(12i) is the characteristic function of the 


set Di, where 12i is a subset of fl, consisting of points (xi,X 2 ) such that —a/Sxi < 
X 2 < y/ixi. Moreover, $ satisfies the boundary conditions (l20ll . Fig.l presents the 
velocity field for this case. The function A was chosen such that the pressure ttq 
in the center is lower that on the boundary of 12, therefore the vorticity is cyclonic 
(anticlockwise). A couple of lines of discontinuities inside the vortical structure can 
be used for modeling of cold and warm fronts in a cyclone. 


5. Estimate of the discrepancy 


5.1. Smooth case. 

Proposition 5.1. Let the stream function <i>, the solution to |27| ) and (031, be 
differentiable in Then the discrepancy Q (see m) is completely defined by 

initial data for the bearing field tti . 






FROZEN AND ALMOST FROZEN STRUCTURES 


7 



Figure 1. 

Proof. Since \Q\^ < 4cq sup |V7ri(t,x)p, it is enough to prove that 

xeR2 

sup |V7ri(t,x)p < sup |V7ri(0,x)p. (22) 

xeR2 xeR^ 

Due to the divergence-free condition equation o has the form 

9*771 -I- (Vtti • u ) = 0. 

We take the gradient of this equation and then multiply by Vtti to obtain the 
transport equation with smooth coefficients (ui,U 2 ) 

9t|VTi|2 + (V|V7ri|2.u) =0. (23) 

Thus, the initial value are transported along smooth characteristic curves. This 
immediately implies (|2^ . □ 

5.2. Non-smooth case. Many papers are devoted to the transport equation with 
non-smooth coefficients. R.J. DiPerna and P.-L. Lions have proved this uniqueness 
result under the assumption the coefficients are in the Sobolev class (locally 

in space), and later L. Ambrosio extended this result to coefficients of class BV 
(locally in space). Previous results on two-dimensional transport equation are due 
to Bouchut and Desvillettes [13], Hauray [14], Colombini and Lerner m- 

To show the maximum principle for (I23p in the case on non-smooth velocity 
u, we use the following theorem on the well-posedness of the transport equation 
proved in US], m- 

Theorem 5.1. (|16jl Let b : be a bounded, divergence-free, autonomous 

vector field on the plane admitting a Lipschitz compactly supported potential f : 
R^ —>■ R, that is b = V^f. The Cauchy problem for 

ut -\- div(bu) = 0 

admits a unigue bounded generalized solution in the sense of distributions for every 
bounded initial datum if and only if the potential f satisfies the weak Sard property. 
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Corollary 5.1. Let the stream-function $, taking part of the generalized solution 
to system m, (w be Lipschitz, compactly supported and |V$| ^0 almost every¬ 
where. Then 

||V7ri(t,x)||i=o < ||V7ri(0,x)||z,«., 

therefore the discrepancy is completely defined by initial data for tti . 

Proof. li is enough to write in the divergent form 

9t|V7ri|2 +div(|V7ri|2v_L$) = 0. 

and apply Theorem 15.11 As noticed in [16) . the weak Sard property is implied by 
|V<i)| 0 almost everywhere. □ 
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